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Abstract 

We study the interference of resonant Higgs boson exchange in neutrahno 
production in /i annihilation with longitudinahy polarized beams. We 
use the energy distribution of the decay lepton in the process Xj ~^ to 
determine the polarization of the neutralinos. In the CP conserving Minimal 
Supersymmetric Standard Model a non-vanishing asymmetry in the lepton 
energy spectrum is caused by the interference of Higgs boson exchange chan- 
nels with different CP eigenvalues. The contribution of this interference is 
large if the heavy neutral bosons H and A are nearly degenerate. We show 
that the asymmetry can be used to determine the couplings of the neutral 
Higgs bosons to the neutralinos. In particular, the asymmetry allows to de- 
termine the relative phase of the couplings. We find large asymmetries and 
cross sections for a set of reference scenarios with nearly degenerate neutral 
Higgs bosons. 



1 Introduction 



At a muon collider, neutral Higgs bosons are produced as s-channel resonances 
in iJ,'^n~ annihilation 12 El- Therefore a muon collider is an excellent tool to 
study the properties of a heavy scalar or pseudoscalar neutral Higgs boson. The 
CP conserving Minimal Supersymmetric Standard Model (MSSM) contains three 
neutral Higgs bosons, a light scalar h, a heavier scalar H, and a pseudoscalar A. 
Since the two heavier neutral Higgs bosons may decay into neutralinos a muon 
collider opens the possibility to test supersymmetry through the interaction of the 
Higgs sector with the neutralino sector. 

Neutralinos, the supersymmetric partners of the neutral Higgs and gauge bosons, 
are expected to be light in large regions of the supersymmetric parameter space. A 
scan of their production line shape may allow to separate the contributions from the 
different resonances [3] and a determination of the coupling mechanism of the Higgs 
sector to the higgsino and gaugino sectors j3]. However, the two heavier resonances 
are nearly degenerate over much of the (tan (3, m^) parameter space and a large 
overlap will make the determination of the resonance parameters a difficult task. 
The production process dependence on the polarizations of the initial muons and 
that of the neutralinos may allow to disentangle the contributions of the different 
resonances. 

In this paper we study the 'XaX'j production process with longitudinal 

muon beam polarization, where i and j label the mass eigenstates of the neutralinos. 
The interference of the Higgs boson exchange channels with different CP quantum 
numbers can be sizable when the mass difference of these Higgs bosons is at most of 
the same order of their decay widths [Sj. The neutralino polarization, averaged over 
the neutralino production angles, results from the interference of the neutral Higgs 
bosons since the contributions from Z and fl exchange vanish, due to the Majorana 
character of the neutralinos [H]. In the center of mass system (CMS) the energy 
distribution of the decay lepton of the decay Xj depends on the longitudinal 

neutralino polarization. The interference effect is proportional to the sum of the 
polarizations of the initial fermions, as well as that of the neutralinos, and thus 
vanishes for unpolarized beams. 

This paper is organized as follows: in Section 2 we present the formalism and dis- 
cuss the lepton energy distribution asymmetry. In Section 3 we study the possibility 
of using the asymmetry to determine the Higgs neutralino couplings. In Section 4 
we present numerical results and in Section 5 we present a short summary and draw 
the conclusions. 

2 Definitions and Formalism 

We study in the MSSM neutralino production 



with longitudinally polarized beams for center of mass energies around the reso- 
nances of the heavy neutral Higgs bosons H and A, and the subsequent leptonic 



/i /i 
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two-body decay of one of the neutralinos 

^o^£±£T £ = e,/i,r. (2) 

The production proceeds via the resonant exchange of H and A in the s-channel, 
as well as via the non-resonant exchange of the Z boson and of the light Higgs boson 
h in the s-channel and of the t- and u-channel exchange of fjiL,R- 



2.1 Lagrangian and couplings 

The interaction Lagrangians for neutralino production via Higgs exchange are (in 
our notation we follow closely [1113111) 



2^' Xi [Cij -TL + ^RjXj (p, 



(3) 
(4) 



where Pr,l = |(1 i 7^); fi' is the weak coupling constant and (j) = H, A, h. In the 
neutralino basis {7, Z, H^, ifg } the muon and neutralino couplings to H and A are 

m 



cos a 
2mw cos l3 ' 



tan /?, 



S'', sin a, 



2mw 

-Qij cos a , 

-z(g:;,sin/?-s:;cos/3) 



2 cos 9w 
1 



(5) 
(6) 

(7) 
(8) 

(9) 

2cos^H'^'-"'-^'^''"''^' 
where a is the Higgs mixing angle, tan/5 = V2/V1 is the ratio of the vacuum ex- 
pectation values of the two neutral Higgs fields, 9w is the weak mixing angle and 
is the unitary 4x4 matrix which diagonalizes the neutralino mass matrix Y. If 
CP is conserved Y is real and the matrix can be chosen real and orthogonal: 
NiaYa/sNji^ = riimy,^6ik, where rn^^.i = 1,...,4 are the masses of the neutralinos 
and ?7j = ±1 is related to the CP eigenvalue of the neutralino x^. The muon and 
neutralino couplings to the lighter Higgs boson h are obtained substituting a by 
a + 7r/2 in eqs. © and (0). 

The interaction Lagrangian for neutralino decay into a lepton and a slepton of 
the first two generations is given by 

gfiiPRxfh + 9f!,iPLX% + h.c. 



with couplings 

fej 



Lcost^vi/ 
-V2ei sin 6w 



{Tu - ei sin^ Ow)Nj2 + ei sin OwNji 



tanOwN. 



'i2 



(12) 
(13) 
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where and T^i denote the electric charge and third component of the weak isospin 
of the lepton 

Mixing can safely be neglected for the scalar leptons of the first two generations, 
£ = e,fl. For the neutralino decay into staus Xi ~^ ^nT, we take stau mixing into 
account and write for the Lagrangian ^Oj 

jC-rfx. = grnHal^PR + bl^PL)x- + h.c. , n = 1, 2; J = 1, . . . , 4, (14) 
where the coefficients a^^ and are given in Appendix |X1 

2.2 Cross section and lepton energy distribution 

To calculate the cross section for the combined process of neutralino production 
and decay we use the spin density matrix formalism of ^1], as e.g. for neutralino 
production in e~^e~ annihilation in [Hj. 

The spin density matrices of production and of x° decay are given 

by 

Pxj^j = '^'^\\,Tx*\'^, (15) 

A, 

Px'^x, = T5;tZ. (16) 

where 7^^^ ^"^^ '^x denote the helicity amplitudes for production and decay, re- 
spectively. The amplitude squared for production and decay is then 

\T\' = |A(x;)PEpS.;P?..., (17) 
A, a; 

with the propagator A(x°) = i/[p'^. ^^Xj '^'^^Xj'^Xj]- Here p^., m^^ and V^. denote 
the four- momentum squared, mass and width of Xj', respectively. For the propagator 
we use the narrow width approximation. 

Introducing a suitable set of spin vectors [TT], [HI, given in Appendix IB. 11 the 
spin density matrices, eqs. (|T5|l and (|T6|) . can be expanded in terms of the Pauli 
matrices t"", a = 1, 2, 3, 

3 

Pf^A; = ^A,A^P + $:r;,,E«p, (18) 

a=l 
3 

P?.A, = 5a;a,/^ + E^aV.S^- (19) 

a=l 

In eq. p8p Sp /P is the longitudinal polarization of the neutralino, Sp / P the trans- 
verse polarization in the production plane and Sp/-P that perpendicular to the 
production plane. Inserting the density matrices, eqs. (fT^ and (jlHI), into eq. (fTTjl 
we obtain 

3 

|T|2 = 2\A{x])?{PD + Y.^1,T.l). (20) 

a=l 
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The first term in eq. (j^Uj) is independent of tlie neutralino polarization and tlie 
second term describes tlie spin correlation between production and decay of the 
neutralino x^. 

The resonant contributions from the s-channel exchange of H and A to P and Ep 
are denoted by Pr and S^, respectively. Explicit expressions for Pr and S|. are given 
in the Appendix, eqs. ()B.4jl and ()B.5jl . while the resonant contributions T,]^ and to 
the transverse polarization of the neutralino vanish. The non-resonant contributions 
to the production density matrix coefficients from Z and slepton exchange can be 
found in jUj. The contributions from exchange of the lighter Higgs boson h are 
numerically negligible. The interference of the chirality violating Higgs exchange 
amplitudes with chirality conserving Z and slepton exchange amplitudes is of order 
m^/y/s and can safely be neglected. The coefficients D and E|, for the two-body 
decay of the neutralino into a lepton and a slepton are given in Appendix IB.3t 
eqs. (|R6|) and (|R8|l . for £ = e, /i and eqs. (p^ and (IBTT|1 for i = t. 

The kinematical limits of the lepton energy in the CMS are 



^rnaxC^n) ^ ± A,, (21) 



with 



^ 2 2ml. 



Xj 



2 2ml 

Xj 



\PJ- (23) 



It is useful to introduce the average 

P = Pd^x^^ (24) 

S| = i- y T?pdn^o, (25) 

over the neutralino production angles in the CMS. Then, the integrated cross section 
for neutralino production, eq. ([T)), and subsequent leptonic decay ^^^n; with 

n = R, L for i = e, fi (n = 1, 2 for £ = r), is given by 

1 (i^l - ^V) - 

< = aA 2 2 3 (26) 

Xj -^3 



Explicit expressions for D are given in eqs. ()B.6|) and ()B.9|) for i = e, fi and i = t, 
respectively, and Xij is the triangle function, defined in Section El The energy 
distribution of the lepton is 



dEf 2A, 



Sp (E^ — El 



P 



(27) 



Here r]i± = =f1. Further, 77^^ = +1 and r;^^ = —1 for the decay into cr^JIr and 
eL,P^L, respectively. For the decay x° ^"^^1,2 ^^e lepton energy dependent term 
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50 100 150 200 

Ei[GeV] 

Figure 1: Normalized primary and secondary lepton energy distribution, with Af = 0.2, 
m^o = 250 GeV, m,- = 200 GeV, m^o = 60 GeV and Js = 450 GeV. The dashed curves 

correspond to the decay chain X2 ~^ ^~^rj ~^ Xi^"*" ^^'^ t^is dash-dotted curves to 



in (j2Zj) is suppressed, due to stau mixing, by 





2 _ 




to 




2 + 


<3 


2 



(28) 



where the X^'^n.T couphngs a'^j and b'^j are defined in eq. ()A.2jl . 

Due to the Majorana character of the neutrahnos, the contribution to Sp from 
the non-Higgs channels is forward-backward antisymmetric jO], whereas that from 
Higgs exchange is isotropic. Then, the non-resonant contribution in eq. ()25|) vanishes 
and, neglecting the interference of the resonant amplitudes with the Z and slepton 
exchange amplitudes, 

s| = (29) 

From eq. ()B.5|1 follows that and thus the energy dependent term in eq. ()27|1 . are 
proportional to the interference of the H and A exchange amplitudes. 



2.3 Lepton energy distribution asymmetry 

For the processes — ^ 'XaX^ with subsequent decay x° — > C^C'^i, with £ = e,/i, 
and Xj cis well as the charge conjugated decays, we define the asymmetries 

and A^- , with n = i?, L for £ = e and £ = /i, and n = 1, 2 for £ = r, 

a^±{Ee > Ee) + a]^±{Ee < Ee) ^ ' 

= IveVi^^ (31) 

in order to isolate the term dependent on the interference of H and A in eq. ()27|) . 

The slepton decays subsequently into a neutralino and a secondary lepton. The 
latter needs to be identified from the primary lepton. Therefore, it is useful to define 
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the asymmetry 



equivalent to A^^ since for the primary lepton A^- = —A^+. The advantage of the 
new asymmetry, eq. (jH^ . is that the largest part of the non- irreducible background 
from the secondary lepton drops out because its energy distribution is only weakly 
dependent on the sign of the lepton charge. In fig. ^ we show the normalized energy 
distributions of the primary and secondary leptons for both charge cases, for a 
sample scenario with Af = 0.2. 

Denoting by aR^n^ jjT xix^) the resonant contribution to the production cross 
section a{fi~^fj,~ — > Xi'x?) we relate P to Pr by 



and express A^ in the form 



~ 2^^ a(^V"- X?X?) ^ 
Vf = ^. (35) 

The contribution of the H-A interference to the asymmetry, eq. ()34|) . is contained 
in the coefficient Vj^, which has the following dependence on the longitudinal beam 
polarizations P^ of fi~^ and P^ of 

= ^Ip^pA r-^ (36) 

where Vf^^R^ = Vf for = P^ = 1, i.e. for right handed /i"*" and /i~ beams. 

Since is proportional to the interference of the H and A exchange amplitudes 
a non-vanishing asymmetry of the lepton energy distribution is a clear indication of 
nearly degenerate scalar resonances with opposite CP quantum numbers. 

The statistical significance of the asymmetry A^ can be defined by 



57 = |^^|y'2a(/i+/i--.x°X°)5^(x°-^^"^+)^e//, (37) 

where Ceff = e"-^ denotes the effective integrated luminosity, with e" the detection 
efficiency of the leptons in the processes ^^^n ^^"^ ^ integrated luminosity. 



3 Determination of the Higgs-neutralino couplings 

A measurement of the asymmetry of the primary lepton energy distribution opens 
the possibility to determine the if-neutralino and A-neutralino couplings using 
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eqs. dHU) and (jHH)- Inserting in eq. (fH^ the expressions of Pr and E|,, eqs. ()B.4|) 
and ()B.5|) . we obtain, for right handed /i"*" and fi~ beams (P^ = = 1): 



^ _ 27.,Re(A(i7)A*(A))^4s-. 



where 



A(0) = z[(s-mj) +zm^r<^]-\ </> = i/, A (39) 
4 = ^ - iVim^^ ± Vjm^,? , (40) 

"^'^ = |c(^)c(Am)|' ^^^^ 
^ Im(cSf^cSf-)lm(c(^^)c(^^)-) 

^ij I I I 

Since we assume CP conservation takes the values ±1 for interfering amphtudes 
with opposing CP eigenvalues, as is here the case, and vanishes for interfering 
amplitudes of same CP. Notice that the functions sfj depend on the relative CP 
phase factor of the neutralinos rjij = rjirij, with rji, rjj defined in Section EUl 

We obtain Vj^ from using eq. (jHH). The neutralino-slepton couplings needed 
to evaluate rj^ will have been precisely studied at a linear collider, see, e.g., [IBIIII!, 
and the resonant cross section of neutralino production aR{fi~^fi~ — >• Xi'Xj) can be 
obtained subtracting the continuum contributions from the integrated production 
cross section a{fi~^fi^ X^xf)- The continuum can be estimated extrapolating the 
production cross sections measured below and above the resonance region [TH] . 

It is possible to solve eq. (j38|) for 'jijrij. To determine •jij and r^j the resonance 
parameters, as well as the peak and off resonance cross sections, need to be precisely 
known. The widths and masses of nearly degenerate Higgs resonances with different 
CP quantum numbers may be determined with transverse beam polarization, which 
enhances or suppresses the Higgs boson cross section depending on the Higgs CP 
quantum numbers |16j . 

The product of couplings 

hj = Im(c;f ^cSf *)Im(c(^^)c(^^>) (43) 
can be determined with a measurement of A"cr(/i+/i~ XAX'j)- 

^^a(/xV-X°x;) = ^^.V«(/iV-X°X?)^ 



R 



9 ^J±(^pL J tdL 

where in the last equality we used the relation 



" ^'^^mis^^' + Pt)Re{{A{H)){A{A)r}r^,r^-,h,, (44) 



a^(/iV-X°x;) = ^Pr. (45) 
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Scenarios 


SPSla 


B5 


BIO 


B20 


B5' 


BIO' 


B5" 


BIO" 


tan (3 


10 


5 


10 


20 


5 


10 


5 


10 


M2[GeV] 


192.7 


280 


280 


280 


280 


280 


280 


280 


/x[GeV] 


352.4 


250 


250 


250 


250 


250 


250 


250 


m^o [Ge Vl 


359 


255 


257 


258 


255 


257 


255 


257 


m^o [Ge V] 


177 


209 


212 


214 


209 


212 


209 


212 


m^o [Ge V] 


96 


128 


131 


132 


128 


131 


128 


131 


mo[GeV] 


100 


100 


100 


100 


100 


100 


100 


100 


me^[GeV] 


143 


173 


173 


173 


173 


173 


173 


173 


m^[GeV] 

L J 


393.6 


450 


450 


450 


350 


350 


550 


550 


mjy[GeV] 


394.1 


451.4 


450.4 


450.1 


351.9 


350.5 


551.1 


550.3 


rA[GeV] 


1.38 


2.29 


2.11 


4.33 


0.43 


0.82 


3.63 


3.34 




0.93 


1.12 


1.33 


3.68 


0.27 


0.71 


2.83 


2.76 



Table 1: Reference scenarios. GUT relations are assumed for the gaugino soft 
breaking mass parameters Mi = 5/3 tan^ 9wM2 and for the slepton mass parameters 
|17j . The resonance parameters are evaluated with HDECAY |18j . 

the triangle function is defined by Ajj = and Pr and are given in eqs. ()B.4|) 

and ()B.5|1 . respectively. 

Notice that the H-A interference in neutralino production depends on 7jj, while 
pure H oi A exchange does not. Therefore, a measurement of the lepton energy 
asymmetries provides unique information on the Higgs-neutralino couplings. 

4 Numerical results 

We present numerical results for the neutralino production cross sections a(/i^/i^ — > 
the asymmetries Af and A}^ of the lepton energy distribution and the sta- 
tistical significance at center of mass energies around the resonances of the neutral 
Higgs bosons H and A. We study the dependence on the MSSM parameters tan /3, 
/i, M2 and mA in the mixed scenarios B and in the gaugino-like scenario SPSla 
defined in Tabled Further, we discuss in scenario SPSla ^2] the influence of beam 
polarization. 

In order to reduce the number of parameters we assume the GUT relations for 
the gaugino mass parameters. Mi = 5/3 tan^ 6'vyM2, and for the slepton masses of 
the first two generations, i = e, n, ^Zj 

m| cos 2/3 sin^ 6'vi/, (46) 
m|cos2/3(^ - sin^^H-). (47) 



= mg + + O.23M2 - 
m? = ml + m'j + 0.79M| - 
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Scenarios 


SPSla B5 BIO B20 


Ar[GeV] 


-254 


BR(x^^r4)[%] 


3.2 16.3 15.2 11.3 


BR(x^^r-f+)[%] 


42.5 17.3 19.6 27.4 



Table 2: Neutralino branching ratios, £ = e,/i. 

where mz is the mass of the Z boson and mo is the scalar mass parameter. The 
masses of the stau mass eigenstates fi and f2 are obtained diagonalizing the stau 
mass matrix, eq. ()A.5j) . Our scenarios have been chosen such that the mass of the 
lightest stau, fi, is of order m^^, and the mass of the heavier stau, f2, is of order 
m^^, and m^^ < m^_. < rrii^. Therefore, neglecting the three-body decays of only 
— > Hr, £ = e, /i and — > rfi contribute to the energy spectrum of the leptons. 
We show the neutralino branching ratios into lepton slepton pairs in Table |21 Details 
of stau mixing can be found in Appendix El 

4.1 X1X2 production 

We first discuss, for X1X2 production, the dependence of the asymmetries and cross 
sections on the beam polarization in scenario SPSla, the dependence on tan/3 in 
scenarios B5, BIO and B20 and the dependence on rriA in scenarios B5',B5", BIO' 
and BIO". All the scenarios are defined in Tabled 

4.1.1 Beam polarization 




385 390 395 400 385 390 395 400 385 390 395 400 



^[GeV] v^lGeV] ^[GeV] 

Figure 2: ^t+^u" XiX2^ X2 ^ ^~^r^ ^ = e,/i, for scenario SPSla. (a): Af , (b): 
neutralino production cross section and (c): significance with luminosity times detection 
efficiency e£ = Ceff = 0.5/6^^ (for X2 ~^ ■^~^r)^ beam polarizations: = = —0.2 
(dash-dotted), -0.3 (dashed), and -0.4 (solid). 

In figs. |2t and |2b we show the asymmetry Af, i = e, fi, and the cross section 
a{fi^fi~ — >• XiX^)' respectively, for scenario SPSla as a function of the center of 
mass energy around the heavy Higgs resonances. Since the resonances are completely 
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overlapping the interference between the CP even and CP odd amphtudes is large, 
resulting in large asymmetries in the resonance region. The largest asymmetries are 
found at ^/s ~ rriH where the CP even and CP odd amplitudes are of the same order, 
because, due to the relative CP phase factor 7712 = ?7i?72 = 1 of the neutralinos, the 
CP even amplitudes are P-wave suppressed. The largest production cross sections 
are found at ^/s ~ itia- 

We show the cross section, the asymmetry and the significance of scenario SPSla 
for P^ = = —0.2, —0.3, —0.4. The dependence on the longitudinal beam 
polarization of the resonant cross section is given by the factor 1 + P^Pf, and is 
rather weak for the polarization degrees expected at a muon collider. From eqs. (j3l|l 
and follows that the polarization dependence of the asymmetry is then roughly 
Ai ~ P+ + P^. For the statistical significance, defined in eq. (pTfjl . follows then 
57 ~ + Pf. 

4.1.2 Stau mixing dependence 

The asymmetry for the r energy spectrum depends strongly on the mixing in the 
stau sector. The r energy asymmetry is obtained by Al = Tj^Af, eqs. flMj) and 
(j^H|) . For the SPSla scenario 77^ = —0.87. Notice that the asymmetries A^. and 
Af^ have opposite signs. The marked difference between Al and Af is due to stau 
mixing, which allows the lightest scalar tau fi to have a large left component. For 
the gaugino-like SPSla scenario the second lightest neutralino X2 is wino-like, and 
thus has large left handed couplings to lepton-slepton pairs. Therefore, X2 decays 
dominantly into r-fi pairs, see the branching ratios for SPSla in Table |21 For 
Aj. = /itan/3 the stau mass matrix is diagonal, eq. ()A.5|) . Then, the branching 
ratios for the decays X2 ~^ ^^R X2 ~^ ''"^1 comparable in size. However, A}^ 
is still smaller than Af, with 77^ = 0.53, due of the larger couphngs to the higgsino 
components to the scalar taus. 

In the resonance region we find, for P^ = P^ = —0.3, Sf ^ l.b^J Ceff[fh~^], 
£ = e,/i, and S]. ~ A.b^J Cefflfh^"^]. In fig. Efc we show the statistical significance, 
defined in eq. (jSZ|), for Af, i = e,fi with an effective integrated luminosity Ceff = 
0.5fb-\ 

4.1.3 tan/3 dependence 

In fig.lHK we show the asymmetry Af, i = e, fi, for scenarios B5, BIO and B20, for 
P^ = P^ = —0.3. These scenarios differ only by the value of tan/3. For increasing 
tan (3 the mass difference mn — mA decreases and the widths Th and increase. 
This results in a larger overlap of the resonances which leads to large asymmetries 
in the resonant region. For tan/3 = 5, with only partial overlap of the resonances, 
the asymmetry is further suppressed by the relative larger continuum contribution 
to the cross section due to the smaller Higgs-muon couplings. However, it shows 
an interesting energy dependence due to the different complex phases of the Breit- 
Wigner propagators of H and A. The maximum of the asymmetry is found at 
y/s ^ rriH, as already discussed for scenario SPSla. 
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450 455 460 



450 455 460 




450 455 460 



Figure 3: /i"*"//" X2 ~^ ^ = e,/U for scenarios B5, BIO and B20. (a): Af, 

(b): neutralino production cross section and (c): statistical significance witli luminosity 
times detection efficiency eC = C^ff = 0.5/6"^ and P;^ = = —0.3. tan/3 = 5 (dash- 
dotted), 10 (dashed), and 20 (solid). 



In fig.Eb we show the cross sections (T(yU^/i^ X1X2) scenarios B5, BIO and 
B20. The largest peak cross sections are found for tan/5 = 10. For tan/3 = 5 the 
resonant cross sections are suppressed by the smaller Higgs-muon couplings, while 
for tan /3 = 20 they are suppressed by the larger resonance widths. 

In the resonant region we find, for i = e, jj, > S^J Cefflfb^^] for tan/3 = 
10 and tan/3 = 20, while for tan/3 = 5 the statistical significances reach ~ 
^ Ceff[fb~^ at ^/s ~ rriH- In fig. Efc we show the statistical significances for an 
effective integrated luminosity C^ff = 0.5/6~^. 

The effect of stau mixing on the asymmetry A!,. = "rjlAf, with 77^ defined in 
eq. (PHj) . increases with tan/3. It is weaker in the mixed scenarios than in the 
gaugino-like SPSla scenario, as can also be observed comparing the neutralino 
branching ratios of Table |2l We find, for scenarios B5, BIO and B20, rj^ = 0.88, 
0.49 and -0.31, respectively. The statistical significance for A^ is obtained from 
eq. (jSZj), where the branching ratios of X2 ~^ ^ = ^1 l^i ^^nd X2 ~^ Pf^ are 

shown in Table |21 

For At- = /itan/3, i.e. for a diagonal stau mass matrix with fi = f^, we find 
rjl = 0.96, 0.79 and 0.30 for scenarios B5, BIO and B20, respectively. 

4.1.4 rriA dependence 

In fig.E^ we compare the asymmetries Af, i = e, fi, for scenarios B5', B5 and B5", 
with different values of m^, as a function of ^/s—mA, for = = —0.3. In fig.^Jj 
we show the corresponding cross section a{fi^fi~ — XiX2)- -^^^ larger Higgs masses 
their widths increase, and thus the interference of the H and A exchange amplitudes. 
However, the asymmetries are reduced by the larger continuum contribution to the 
cross section. 

For smaller Higgs masses, here for rriA = 350 GeV, threshold effects are stronger. 
Since 7712 = 1, the asymmetries nearly vanish for ^/s ~ rriA, where the largest cross 
sections are found, while the largest asymmetries are found at ^/s ~ rriH- The 
asymmetries change sign between the two resonances, due to the complex phases 



12 



0.2 
0.15 

0.1 
0.05 



-0.05 
-0.1 

-0.15 
-0.2 



1 1 

-(a) 


1 1 




1 1 



-4 -2 2 4 6 
^-mA[GeV\ 



1000 




-4 -2 2 4 6 



A': 




2 4 



1000 




2 4 



Figure 4: X1X21 X2 ~^ ^ = ^'A*' ^^^^ = = -0.3, for scenarios B5, 

BIO with ruA = 450 GeV (dashed), B5', BIO' with niA = 350 GeV (sohd) and B5", 
BIO" with niA = 550 GeV (dotted), (a): Af' for tan/? = 5, (b): neutrahno production 
cross section for tan/3 = 5, (c): Af' for tan/3 = 10, (d): neutrahno production cross section 
for tan /J = 10. 



of the propagators, and the maxima of \Af\ are found at center of mass energies 
slightly above and below rriH and not on top of the CP even resonance. For larger 
values of niA, here for = 550 GeV, the peak cross sections are suppressed by the 
larger widths. To a lesser degree, they are enhanced by the larger phase space for 
neutrahno production. 

In figs. |3Jd and|3Ji we show the analogous figures for scenarios BIO', BIO and 
BIO" , with tan (3 = 10. The effect of larger Higgs masses is weaker than for tan /3 = 
5 because the overlap of the resonances is already large for rriA = 450 GeV. 

The statistical significances at the center of mass energies where \Af'\ is maximal 
is ~ O.S^y Ceff[fb~^] for scenario B5' and ~ lAy^ Ceff[fb~^ at ^/s ~ rriH for 
scenario BIO'. 



4.2 XiX2.X2X2 and xiXa production 

In fig. El we show, for scenario BIO, the asymmetry Afji = e,/i, and the cross sec- 
tions a{fi^fi~ — > Xi'Xj*) for the kinematically allowed pairs X1X25 X1X3 X2X2' 

= P^ = —0.3. Threshold effects are stronger in X2X2 X1X3 production than 
in X1X2 production. Therefore the P-wave suppression of the CP even (CP odd) 
amplitude for X2X25 iXiXs) production is stronger, since 7722 = 1 (^713 = —1)- The 
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Figure 5: XiXp Xj ^ ^ ^w^ — ^' scenario BIO with i = I, j = 2 (solid), 

i = 1, j = 3 (dash-dotted) and i = 2, j = 2 (dashed), (a): Af and (b): neutrahno 
production cross section, with Pf = = —0.3. 



asymmetries for X1X2 X2X2 production, however, are comparable in size, since 
the continuum contribution to X2X2 production is very small. For X1X3 production, 
the asymmetry is smaller due to the interplay of the widths and the amplitudes, 
with Th < Fa and A exchange suppressed, which results in a smaller interference of 
the two amplitudes. Notice also the different energy dependence of the asymmetry, 
with maxima at ^/s < rriA and ^/s > mn, and of the cross section, with a maximum 
at y/s ~ niH- 

For the asymmetry in X2X2 production the statistical significance for i = e, fi 
in the resonance region is ~ 2^^ Ceff[fb~^]. For X1X31 the statistical signifi- 
cance for £ = e,yU. is significantly smaller, of order Sf^ ~ 0.6^£e//[/^~^] in the 
resonance region, because the branching ratios of xt iiito lepton and slepton pairs 
are strongly suppressed by the competing decay channels X3 — * Zxi and X3 hxi, 
with BRixl iin) = 1% for £ = e, /i, and BR{xl ^ rh) ~ 5%. 

4.3 /i — M2 plane 

The MSSM parameters /i and M2 affect strongly the neutralino couplings both to 
the Higgs bosons as to the lepton-slepton pairs. The lepton energy asymmetries for 
neutralino decays into leptons of the first two families, e and n, do not depend on the 
couplings, while the dependence of stau mixing on the neutralino character has been 
briefly discussed in Section 14.1. HI The neutralino couplings to H and A are both 
enhanced in mixed scenarios since Higgs bosons couple to a higgsino-gaugino pair. 
Therefore, the Higgs boson widths, and thus the interference of the resonances, are 
also enhanced. In figs and IHb we show contours in the /i — M2 plane of constant 
7i2ri2 and Im(c|^^c^^'**), respectively, for tan /3 = 10 and tha = 450 GeV. Notice that 
712 is negative in most of the experimentally allowed parameter space. Therefore 
the sign of the asymmetries for the first two lepton families constitutes a test of the 
Higgs-neutralino couplings in the MSSM. 

The same qualitative dependence of 7i2'"i2 and ku on fi and M2 is found for 
different values of tan/3 and m^. 
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Figure 6: Contours of constant 7i2?"i2 (a) and Im(c.^^c.^^*) (b), for tan/3 = 10 and 

TUA = 450 GeV. In (b) we show lm{c\fc\f*) = 0.01 (dashed), 0.03 (dash-dotted) 
and 0.05 (dotted). The wiggly hues in both figures indicate the level crossing of the 
X2 and Xs states, with rju = 1 in the area below and 7712 = —1 in the area above the 
level crossing line. The gray area is experimentally excluded by m^± < 103 GeV. 

5 Summary and conclusion 

We have discussed the interference of the CP even and the CP odd amplitudes of 
neutral Higgs boson s-channel exchange in fi~^fi~ XiX^ with longitudinally po- 
larized beams in the CP conserving MSSM. To study this interference we use the 
energy distribution of the lepton from the decay i'^i^, £ = e, fi, r. The asym- 

metry of the lepton energy distribution depends on the longitudinal polarization of 
the neutralinos, averaged over their production angles. Since the neutralino longi- 
tudinal polarization is correlated to the longitudinal beam polarization when the H 
and A exchange amplitudes interfere, and vanishes otherwise, the asymmetries can 
be used to determine the couplings of the H and A bosons to the produced neu- 
tralinos times, respectively, their couplings to the muons. In particular, the sign of 
the asymmetry is sensitive to the sign of the product of couplings of the neutralinos 
and of the muons to the Higgs bosons. 

For a set of scenarios we have analyzed the lepton energy asymmetries for Xj 
production with subsequent two-body decay, with emphasis on X1X2 production. We 
find large asymmetries for nearly degenerate heavy neutral Higgs bosons and inter- 
mediate values of tan/3 and niA- Especially for X1X2 production we find statistical 
significances which would allow to measure the asymmetries at a muon collider. 

6 Acknowledgments 

We thank O. Kittel for valuable comments and discussions. This work was supported 
by the 'Deutsche Forschungsgemeineschaft' (DFG) under contract FR 1064/5-2. 



15 



Appendix 



A Stau-neutralino couplings 

The stau-neutralino couplings, defined by the interaction Lagrangian 

Crfx, = gfnfialjPR + hljPL)x] + h-c. , n = 1, 2; J = 1, . . . , 4, (A.l) 

are [H]] 

with TZnm the stau mixing matrix defined in eq. (jA.7|l and 

In eq. (IA.3jl . and //^- are defined by eqs. (fT^ and (|T!?jl . respectively, and 

h^. = {hf^y =m,/{V2mw cos f])N;s, (A.4) 

with rriw the mass of the W boson, rrir the mass of the r-lepton and the neutralino 
mixing matrix in the 7, Z, H^, H2 basis. The masses and couplings of the r-sleptons 
follow from the — tr mass matrix: 

/ m?^ -nirAr \ / ^i? \ 

y -nirAr m\ J \Tl J 

with m?^ and m?^ given by eqs. PU)) and (jTrf) replacing by m^, and 

A^ = A^-;utan/3, (A.6) 

where A^- is the trilinear scalar coupling parameter. The f mass eigenstates are 
(ti7 T2) = {tr, Tl)'TV' with the stau mixing matrix 

cos Of sin Of 
= I I , (A.7) 

- sin Of cos 6*^ 



and 

m? — m| ^ m^A^ , . , 

cosgf = ^ , sin^^ = —= " " . (A.8) 

The mass eigenvalues are 

= ^ (("^'l + "^L) T y/(m?^ -m?^)2 + 4m2A2) . (A.9) 
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B Density matrix 



B.l Polarization vectors 

The polarization vectors s^(Xj) (c^ = l!2,3) of the neutrahno Xj in the CMS are 
chosen so that s^. and s^. are perpendicular to the momentum of the neutralino p^^ 
and s^. is parallel to p^^ . In the reference frame where the four-momentum of the 
neutralino is given by 

p^^, = (^^^-0,0,1^^^), (B.l) 
and where the normal to the production plane is 

P'^-^'P^^ -(0,1,0), (B.2) 



we define the polarization vectors 

4^ = (0; 1, 0, 0), 4^ = (0; 0, 1, 0), sf^ = ^(bl,J; 0, 0, E^^). (B.3) 

B.2 Production density matrix: H and A exchange 

The interaction Lagrangians are given in eqs. Q and (j3)). The contribution of the 
CP-even H and CP-odd A exchange channels to the production density matrix 
coefficients P and Sp, defined by eq. (fT8|l . are 

Pr = i(2-<5.,)(l + P^Pf) 

g^ll^flY^^)l^lA{H)\^ss^^ (B.4) 

S| = (2-5,,)(^^ + ^-) 

^?^Im(c;f)c;f*)Im(c(^^)c(^'^)*)Re{(A(if))(A(A))*}sv/A;;r/,. (B.5) 

Here A{H), A{A) are the Breit-Wigner propagators defined in eq. (j39|l . the functions 
sfj are defined in eq. pUjl and Xij = sfjS~j is the triangle function. 

B.3 Neutralino decay density matrix 

The interaction Lagrangians are given in eqs. (fTTjl and p4j) . The neutralino decay 
density matrix, eq. (HH), is 

a 

The expansion coefficients D and for the two-body decay into a positive charged 
lepton of the first two generations and a right or left slepton are [T^ : 

D = yl/^fK-"^|), (B.6) 
= v79Vi,\'^xA-'l,-Pi)Vj, n = R,L, (B.7) 
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respectively, with s'^. the neutrahno spin- vector defined in Section lB.il pi the lepton 
four- momentum, and r]^ = 1 and r]^ = —1. The couphngs f^j are defined in eqs. fT^ 
and (Uni). 

In the CMS 



= -^.V|/;;r^U- ' E^^jrj,, n = R,L. (B.8) 

For the decay x'j ~^ '^^^n i n = 1,2 the coefficients are 

D = y(|a:/ + |&:/)K-m^), (B.9) 



The couphngs a^^ ,n = 1,2 are defined in eq. ()A.2|) . 
In the CMS 



2 / 2 2 

. / mt . — mi 



= 9\K\' - K\')j^^ [E, - ^^-^^ - E^^j V,, n=l,2. (B.li: 

The decay density matrix coefficients for the charge conjugated processes, Xj 
£~i^,n = R,L an d T~f^,n = 1,2 are obtained inverting the sign of in 

eqs. jHIj), (jEHI), (jBlnll and (pjT]l 
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